Electron dynamics inside a vacuum tube diode through linear differential equations 
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In this paper we analyze the motion of charged particles in a vacuum tube diode by solving linear differential 
equations. Our analysis is based on expressing the volume charge density as a function of the current density 
and coordinates only, i.e. p = p( J, z), while in the usual scheme the volume charge density is expressed as 
a function of the current density and electrostatic potential, i.e. p = p{J,V). Our approach gives the well 
known behavior of the classical current density proportional to the three-halves power of the bias potential and 
inversely proportional to the square of the gap distance between the electrodes, and does not require the solution 
of the nonlinear differential equation normally associated with the Child-Langmuir formulation. 



I. INTRODUCTION 



The motion of charged particles accelerated across a gap is of wide interest in fields such as high power diodes and vacuum 
microelectronics. Child and Langmuir first studied the space charge limited emission for two infinite parallel plane electrodes 
at fixed voltage Vq in vacuum separated by a distance I^.EI^The charges produced at the cathode are made to increase (e.g. by 
increasing the temperature in a thermionic diode or by increasing the power of a laser in case of a photocathode). Interestingly, 
a saturation is observed and the current is said to be space charge limited. A useful approximation for the amount of current 
flowing in such cases is the Child-Langmuir expression for transmitted current. For electrodes having a potential difference V 
and separated by a distance D, the Child-Langmuir law is obtained by solving Poisson's equation 



P_ 
eo 



(1) 



where V is the electrostatic potential, p is the volume charge density and is the permittivity of free space.l^Because this diode 
is assumed to have infinite extension in the x-y directions, we can define the cuiTent density by 



J(z) = p{z)v{z) = -JcL 
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where v is the velocity of the electrons. By charge conservation the current density can not vary with z, hence the current density 
is constant. Now we can find the velocity of the electrons by conservation of energy 
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where m and e are the electron's mass and charge, respectively. In Eq.([3]) we have assumed that the electron is initially at rest 
in the grounded cathode. Solving Eq.(|3]l for the velocity and substituting in Eq.Q we obtain the volume charge density as a 
function of the cuiTent density and the electrostatic potential 
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Substituting Eq.(|4]) into Eq.([T]i we have a second order nonlinear differential equation for the electrostatic potential 
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with the following boundary conditions 



Using the following ansatz for Eq.(|5]l 
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and substituting Eq.(|7]l into Eq.Q we find that the non linear differential equation is satisfied for 



p = - and JcL 



The volume charge density in the gap is determined by Eq.Q 
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The main result is the Child-Langmuir law which states that the behavior of the current density is proportional to the three-halves 
power of the bias potential and inversely proportional to the square of the gap distance between the electrodes. 
Since the derivation of this fundamental law many important and useful variations on the classical Child-Langmuir law have been 
investigated to acc ount for special geometries,"*^ * relativistic electron energies,"^ non zero initial electron ve loci ties, 'S'Sj quantum 
mechanical effects nonzero electric field at the cathode surface,'^ and using vacuum capacitancel^l 

Given the wide interest and applicability in this topic we present the analysis of the motion of charged particles in a vacuum tube 
diode by solving linear differential equations. Our derivation is based on expressing the volume charge density as a function of 
the current density and coordinates only, i.e. p — p{J, x). Our approach gives the well known behavior of the classical current 
density proportional to the three-halves power of the bias potential and inversely proportional to the square of the gap distance 
between the electrodes, and does not require the solution of the nonlinear differential equation normally associated with the 
Child-Langmuir formulation. 



II. ALTERNATIVE APPROACH 



In the Child-Langmuir derivation we have expressed the volume charge density as a function of the current density and the 
electrostatic potential, let us now try to express the volume charge density as a function of the current density and the spatial 
coordinate. The velocity of a charged particle subjected to an electric field is governed by Newton's second law 



TO — = —eE(z) 
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where e — 1.6 x 10 ^^C and E{z) is the electric field. We can obtain the electric field inside the gap by integrating Gauss's law 
from z = to an arbitrary position in the gap 
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where we have used the fact that E{z = 0) = 0. Integrating by parts Eq.( 1 1 1 we have 

E[z) = —z [ z-^dz 
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Let us work now in the adiabatic regime and neglect the contribution of p'{z) in Eq.( 12 1, i.e. the volume charge density is slowly 



varying in space. Substituting Eq.( 12 1 into Eq.( 10 1 we have 
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Integrating Eq.( 13 i from z = to an arbitrary position we have 
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where we have integrated by parts the right hand term. Applying the adiabatic condition to Eq.( 14 1 we obtain 
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Solving for the velocity in Eq.Q and substituting into Eq.( 15 1 we have 
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Using Eq.([T6]l we can express the volume charge density as a function of the current density and coordinates as 



Substituting Eq.([T7| into Eq.([T]) we have 
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Equation ( 18 1 is a linear second order diferential equation which can be solved by quadratures. Integrating Eq.(|T8]l twice and 
using the initial conditions given in Eq.(j6]l we obtain the solution for the electrostatic potential 
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If we apply the boundary condition of the fixed potential anode V{z = 13) = Vo in Eq < 19 1, we can find the current density in 
the adiabatic regime which is given by 
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Equation ( 20 1 is remarkably close to the result for Jc l derived in Eq.([8]), i.e. \ J\ w OAlJc l ■ Substituting Eq.( 20 1 into equations 
([T7]iand(|l9f we find 
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Equations (21 1 are the same exact results derived with the Child-Langmuir formulation. 

If we want to go beyond the adiabatic approximation, we integrate again by parts Eq.([T2| to obtain 
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Assuming now that we can neglect the last term in Eq.(22 1 we have in Newton's second law 



Integrating by parts Eq.(23 i we have 



Neglecting the last term in Eq.( 24 1 we have 
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Solving for the velocity in Eq.Q and substituting into Eq.p5]) we end up with 
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Equation (26 1 is a linear differential equation over p^. Solving Eq.(26 1 for the charge density we find 
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where C is a constant of integration. If we apply the following condition p'{z — D) =0 in Eq.(27 1, i.e. the charge density is 
minimum at the anode, we have the following charge distribution in the vacuum tube diode 



p{z) = 
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Substituting Eq.(28 1 into Eq.([T]) we have 
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Equation ([29]l is a linear differential equation which can be solved by quadratures. Integrating Eq.( 29 1 with Mathematica we find 
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where 2^1 (a, b; c; z) is the Hypergeometric functionP^If we apply the boundary condition of the fixed potential anode V{z 
D) = Vq in Eq.( 30 1, we can find the current density which is given by 
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Evaluating Eq. (31 1 numerically we find that \ J\ w O.hlJcL- Our approach can be easily extended to include the next order 
derivatives in the volume charge density to improve the accuracy of the current charge density, unfortunately this will lead us 
to solve non linear differential equations to express the volume charge density as a function of the current density and spatial 
coordinate and the expressions for the electrostatic potential and charge distribution become far more complex. 
In Fig. ([T]) we have plotted the electrostatic potential in the adiabatic regime and we have plotted the difference between the 
electrostatic potential in the adiabatic and non adiabatic regime, respectively. Note how there is no difference at z = and 
z = Dior the electrostatic potential function in order for the applied bias to be Vq in the vacuum tube diode. In Fig. (jUj) we have 
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FIG. 1: Plots showing (a) the electrostatic potential in the adiabatic approximation with D = 1, (b) the difference between the electrostatic 
potentials in the adiabatic and non adiabatic approximation. Note how the functions for the electrostatic potential have the same value at 2 = 
and z = D such that the applied bias in the vacuum tube diode is Vb and are really close to each other for < 2: < D. 



plotted the volume charge density for the adiabatic and non adiabatic regime where we can see that the functions are practically 
identical except at z — D where the volume charge density for the non adiabatic regime has an horizontal slope. 
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FIG. 2: We have plotted the graphs of the volume charge density for the adiabatic and non adiabatic regime with D — 1. Note how the 
functions superimposed except at z = D where the volume charge density for the non adiabatic regime has an horizontal slope. 
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III. CONCLUSIONS 

The alternative method presented in this article of deriving the V'^/'^ / D"^ scaling of the Child-Langmuir law avoids the need 
of solving a nonlinear differential equation and presents a new insight into the way of approaching the problem of the charge 
dynamics inside a vacuum tube diode. In the Child-Langmuir formulation, the Poisson equation is solved with a volume charge 
distribution that involves the electrostatic potential in contrast with our procedure which provides the volume charge distribu- 
tion as a function of coordinates, so that the solution can be readily determined by simply integrating Poisson's equation and 
applying the boundary conditions at the cathode and anode, respectively. In addition, we have obtained solutions for the elec- 
trostatic potential and volume charge density in the adiabatic and non adiabatic approximations. We have shown that one has 
the same solution for the electrostatic potential and volume charge density as in the Child-Langmuir formulation in the adiabatic 
approximation but with a different current density. 
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